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The influence of a surface in the non-retarded interaction between two atoms 


Reinaldo de Melo e Souza/ W.J.M. Kort-Kamp/ F.S.S. Rosa/ and C. Farina^ 

^Instituto de Fisica, UFRJ, CP 68528, Rio de Janeiro, RJ, 21941-972, Brasil 

In this work we obtain analytical expressions for the non-additivity effects in the dispersive inter¬ 
action between two atoms and perfectly conducting surface of arbitrary shape in the non-retarded 
regime. We show that this three bodies quantum-mechanical problem can be solved by mapping it 
into a two-bodies electrostatic one. We apply the general formulas developed in this paper in several 
examples. Firstly we re-derive the London interaction as a particular case of our formalism. Then we 
treat two atoms in the presence of a plane, re-obtained the result displayed in the literature. After 
we add some new examples. A particularly interesting one is two atoms inside a plate capacitor, a 
situation where non-additivity is very manifest since the plates lead to the exponentially suppression 
of the interaction of the atoms, provided the atoms are separated by a distance of the order of the 
distance between the plates or greater. Our results holds even in the presence of other atoms inside 
the plate capacitor. As a last example we deal with two atoms in the presence of a sphere, both 
grounded and isolated. We show that for realistic experimental parameters the non-additivity may 
be relevant for the force in each atom. 

PACS numbers: 34.35.-|-a,31.30.J-, 34.10.-l-x, 34.20.-b, 34.20.Gj, 31.15.am 


I. INTRODUCTION 

One of the most interesting features of the van der 
Waals dispersive interactions is their non-additivity[l|- 
pointed out for the first time in@ (for a pedagog¬ 
ical exposition see 0 )- It means that the interaction 
between three bodies doesn’t follow from the superpo¬ 
sition principle, or, equivalently, that the presence of a 
body influences the interaction of others. Although such 
effects have been known for many decades, their con¬ 
sequences are not yet fully investigated and there are 
still a lot of research trying to figure them out. One 
of the fertile soils that is being explored is the so-called 
Efimov quantum stateQ, in which a resonant two-body 
force between identical bosons can produce bound states 
in a three-body system even if there isn’t any corre¬ 
sponding two-body bound state. The utilization of Fes- 
hbach resonances to tune the interaction in ultracold 
atom systems has allowed to experimentally probe the 
Efimov states @-[111. These are part of a larger pro¬ 
gramme which studies universal properties of few-body 
systems with large scattering length. [l^-(TB|. It is also 
well-known that the interaction between two Rydberg 
atoms can inhibit all but a single collective Rydberg ex¬ 
citation, a phenomenon called dipole blockade H-iia. 
It has recently been shown that the inclusion of a third 
Rydberg atom can break the dipole blockadefl^ due to 
non-additivity effects. Another very interesting conse¬ 
quence of non-additivity is the great enhancement of the 
dispersive force between atoms in the vicinities of a one¬ 
dimensional transmission line(l@. 

In the literature there are numerous papers dedicated 
to study the non-additivity [l^-[13 ™ some simple sys¬ 
tems. The particular case of the three-atom problem is 
presented in many textbooks 0-0- Unfortunately, due 
to major calculation difficulties, there are not many cases 
analysed, specially those involving macroscopic bodies. 
A method developed by C. Eberlein and R. Zietal[23| 


enables us to evaluate the non-retarded dispersive inter¬ 
action between one atom and a perfectly conducting sur¬ 
face of arbitrary shape, requiring only the knowledge of 
a classical Green function that can be obtained from an 
electrostatic problem. This method has been applied in 
a variety of interesting problems [13 - [13 ■ In this paper 
we generalize this method in order to have two atoms in¬ 
teracting with a surface, and thus we are able to obtain 
analytically the influence exerted by a given perfectly 
conducting surface on the interaction of two atoms in the 
non-retarded regime. As a particular case, we show that 
in the absence of surfaces we recover the well-known in¬ 
teraction between two-atoms, namely, London’s formula. 
This paper is organized as follows. In the next section 
we generalize Eberlein-Zietal’s method to include a sec¬ 
ond atom in the system. We then identify the general 
expression for the non-additivity term and show how it 
is related to an electrostatic problem of a single charge 
in the presence of the conducting body. In section III we 
evaluate the influence of a surface on the interatomic in¬ 
teraction and in the subsequent section we analyse some 
examples. Firstly we re-obtain the interaction energy for 
two atoms in the presence of an infinite conducting plane. 
Then we analyse two atoms inside two parallel infinite 
planes. This example is relevant experimentally[3lj and 
we show that the non-additivity effects are strongly per¬ 
ceivable, leading to an exponentially suppression of the 
van der Waals force between the atoms, provided they 
are kept apart by distances of the order of the separation 
between the plates or larger. This has far reaching con¬ 
sequences. For instance, the van der Waals equation for 
gases presents deviations from ideal gases by taking into 
account the finite volume of atoms and the interatomic 
interactions. Hence, our results show that by placing a 
gas inside a plate capacitor it behaves more ideally. A 
similar exponential attenuation was also obtained in[^ 
in the retarded regime, for atoms inside a rectangular 
box. As a last example, we display the calculation of 


2 


two atoms in the presence of a conducting sphere, both 
grounded and isolated. We leave a final section for con¬ 
clusions and final remarks. 


II. INTERACTION ENERGY FOR TWO 
ATOMS AND A CONDUCTING SURFACE 

In the non-retarded regime the electromagnetic field 
doesn’t have to be quantized. Therefore the interaction 
of an atom and a surface, which is usually dealt within 
a quantum electrodynamics framework, can instead be 
approached by standard quantum mechanics techniques, 
where the interaction hamiltonian to be used in perturba¬ 
tion calculations is given by the instantaneous Coulomb 
interaction!^. The convenience of Eberlein-Zietal pro¬ 
cedure is to relate the quantum mechanical problem to 
an electrostatic one, allowing us to solve the non-retarded 
interaction in the simpler electrostatic domain. 

To begin with, let us consider two atoms A and B 
at positions and rg, respectively, in the presence of 
a grounded surface S. The electrostatic energy of the 
configuration is given by 

[/= iy p(r)$(r)d^r, (1) 

where p(r) gives the charge distribution and $(r) is the 
electrostatic potential, which satisfies Poisson equation 


where Gh satisfies the Laplace equation V^Gh = 0 with 
the boundary condition 


1 

47r|r — r'l 


G//(r,r') 


= 0 . 


res 


( 8 ) 


The equations obeyed by Gh are analogous to those 
satisfied by the potential $i(r) generated by the image 
charges in the electrostatic problem of a charge q at po¬ 
sition r' in the presence of a perfectly conducting surface 
S. If we solve this electrostatic problem we will get Gh 
from the relation 

G^(r,r') = £2^. (9) 


The variable r' is implicitly present in 4>i{v) since the 
image charges depend upon the position r' of the source 
charge. As it will become clear, Gh is the only function 
which must be calculated in order to evaluate the non¬ 
additivity effects of our problem. Therefore, this method 
enables us to effectively replace a quantum mechanical 
problem of two atoms in the presence of a conducting 
body by an electrostatic one of a single charge in the 
presence of the conducting body. To proceed further we 
must specify the charge distribution p{r) appearing in 
eq. ([6]) . We model each atom, in a first approximation, as 
an electric dipole. Hence, the charge distribution is given 
by 


V2$(r) = 

^0 


( 2 ) 


and vanishes at surface S. The solution of equation ([2]) 
can be written in terms of the Green function 


$(r) = — [ G{r,r')p{r')(Pr' 
£o J 


(3) 


where G(r, r') is the solution of the differential equation 
V2G(r,r') = -<5(r-r'), (4) 

subjected to the boundary condition 


G(r,r') 


= 0 . 


(5) 


res 


p{v) = lim g [(5(r - (r^ + h^)) - 5(r - r^)] 

h^->0 

-I- lim g [(5(r - (rs-1-hs)) - ^(r - rs)] 

ghp=dg 

=: pA{v)+pB{r). ( 10 ) 

In the following we write the electrostatic energy of two 
point dipoles in the presence of a conducting surface. 
This allows us to write the quantum hamiltonian interac¬ 
tion for two atoms and a conducting surface by promoting 
d to a quantum operator. Substituting the decomposi¬ 
tion © and eq. (ITUl) into equation ([S]) we obtain 

U = j[pA[^)+PB[r)\G{r,T')[pA{r') +pB[Y')\Sv'd^r 

— Ua Ub H“ Ucrossed i (H) 


By substituting eq-dS]) into eq. o, we write the electro¬ 
static energy as 

U =jcfircPr'p(r)G{r,r')p{r'). (6) 

An immediate particular solution of the equation (|4]) is 
l/47r|r — r'l which, however, doesn’t obey the bound¬ 
ary condition ([S]). This suggests that we decompose our 
Green function in the form 

^ 47r|r-r'| (7) 


where 

/ Pi{r)G{r,r')p,{r')d^Tc'(fr, ( 12 ) 

i = A,B and 

Ucrossed = — f PAir)G{r,r')pBir')d^r'd^r. (13) 

£o J 

To obtain the last equation we used that G(r, r') _= 

G(r', r), whose validity follows from Green’s identity!^. 
Since the Green function is the same in the case of one 
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atom or two atoms, Ui represents the interaction energy 
between point dipole i and the surface S in the absence 
of the other dipole [s^. To unveil the physical meaning 
of Uab let’s employ the decomposition © 


U, 


crossed — 


1 f Pa{v)pb{y') A , ,3 


eo J 47r|r — r'| 




Uab 


+ — [ pAir)GHir,r')pB{r')d^r'(fr.{U) 

eo J 


Uabs 


The first term on the right-hand-side of last equation 
doesn’t depend on the surface. It describes the interac¬ 
tion between two point dipoles in vacuum. The last term 
depends conjointly on both dipoles and the surface. We 
will call it Uabs- Therefore, we write 


U = Ua + Ub + Uab + Uabs ■ (15) 


between the atoms, assumed in the ground state, and 
the surface we proceed perturbatively. In first order of 
perturbation theory we have := {Hint), where (• • •) 
denotes the expectation value of the operator inside the 
brackets in the ground state |0 a,0b) of the atoms. In 
this order only the first two terms in eg. (1181) contribute, 
since for atoms with no permanent dipole moment we 
have 

(0A,0B|dfdf|0A,0B) = {0A\dt\0A){0B\df\0B)=0. 

(19) 

From now on we omit the hats to denote quantum opera¬ 
tors in order to not overburden the notation. Evaluating 
{Ha) and {Hb) we obtain 


ENRiT^A,rB) 




+ 


1 

2so 


Y.{{dg)^)VmVn,GH{vy) 


,( 20 ) 


We see at once one positive aspect of this formalism. 
It enables us to study separately the so-called non¬ 
additivity of dispersive forces, which is totally contained 
in the last term. Substituting eq. m into eq. (HU, per¬ 
forming a Taylor expansion and discarding divergent self¬ 
interaction terms, we obtain Ua and Ub, 


Ua 


Ub 


(d^.V')(dA-V)GB(r,r') 


r—r'—TA 


(ds-V')(dB-V)GB(r,r') 


r—r'—TB 


(16) 


In this way we re-obtain as a particular case Eberlein- 
Zietal formula [ 2 ^ for the interaction of a sole atom with 
a conducting surface. A similar treatment of the terms 
Uab and Uabs gives 


Uab = -(dB-V')(d.4-V) ^ / ,, 

£0 47r|r-r'| 

Uabs = -(ds • V')(dA • V)GB(r,r') 
£0 


r=r^,r'=rB 

.(17) 

r=rA,r'=rB 


Until now, we have been working within classical elec¬ 
trostatics. The passage for the quantum mechanical 
problem is performed by promoting d to an operator 
in eqs. m and (HU, in order to obtain the interaction 
hamiltonian for the quantum mechanical problem of two 
atoms and a conducting surface. Therefore, the previous 
decomposition m can be recast into the form 


Hint = Ha + Hb + Hab + Habs , (18) 


where the operators on the right hand side are obtained 
from UA, Ub, Uab and Uabs by changing d by the quan¬ 
tum mechanical operator d. Note that, as mentioned 
before, the interaction hamiltonian does not involve field 
operators since we are in the non-retarded regime. To ob¬ 
tain the interaction energy for the dispersive interaction 


where we employed {dmdn) = Smn{d‘^), valid for or¬ 
thonormal basis, which are used throughout this paper. 
In other words, in this approximation the atoms don’t 
perceive each other and the interaction of the system is 
the direct superposition of the interaction between each 
atom and the surface S. The non-additivity effects we 
are looking for appear only in the next order. The sec¬ 
ond order contribution to the interaction energy is 


ENRiT^A,rB) 


E' 


I 


{0A,0B\Hint\I){I\Hint\0A,0B) 
Ej - {E^ + E^) 


( 21 ) 

The prime indicates that we must sum over all possible 
states |/) ^ |0yi,0B). Eq and Eq are the energies of 
the ground states of the atoms A and B, respectively. 
Denoting the possible states of atom A by |r) and the 
possible states of atom i? by |s), we may write the above 
formula as 


E^lii^A, Tb) 


E' 


(0.4 , Ob I Hint l^, e) {t ^ s\Hint 1 0.4 , Ob) 
E^ + Eg 


( 22 ) 

where Eg = Eg — Eg and Eg^ = Eg — Eg. Equation 
m shows that in principle we have 16 terms to deal 
with in eq. (1221) but fortunately as we shall see, most of 
them either vanish or are irrelevant for our purposes. For 
the sake of clarity we analyse them separately. To begin 
with, the term 




, {0A,0B\HA\r,s){r,s\HA\0A,0B) 


Egr 




= -E 


, {0A\HA\r){n\HA\0A) 


Egr 


(23) 


since only intermediate states with |s) = |0b) survive in 
the summation. This term depend only on atom A and 
stands for the second-order contribution to the interac¬ 
tion between the atom A and the surface. An analogous 
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term holds for atom B. Since these two terms doesn’t 
contribute to non-additivity we shall henceforth neglect 
them. 

All the other 10 terms involving Hi vanish. Indeed, 
the crossed term involving Ha and Hb, 


= . (29) 

where := — fs- Therefore, substituting last 

equation into (l27ll we get 


-E' 


(0,0|g^|r,s)(r,s|gB|0,0) 

K + Eg 


(24) 


because the term (0,0|iJ^|r, s) isn’t zero only for inter¬ 
mediate states with s = 0, while (r, s|iJB |0, 0) vanishes in 
such cases. From eq. dm) we see that the same argument 
applies to the crossed term 


/ (0,0|iJyi|r,s)(r, s|R^bs|0,0) _ 
^ Eg + Eg - 

1 (0,0|RA|r, s)(r, s|(if dj^|0,0) 

^0^ Eg+Eg 


X 


V,V'GB(r,r') 


r—rA-r'—VB 


(25) 


Summing up, equation (l22ll becomes 


_ {0,0\dfdf\r,s){r,.s\dgd^\0,0)^^ 

^ {4neorR%{Eg + Eg) 

{Sij — S{RAB)i{RAB)j^ (Smn — i{RAB)m{RAB)n^ ■ 

(30) 


As usual, we adopt the notation (0|d^|r) = and, for 
the atom B, (0|d^|s) = (i°®. Therefore we write 

{0,0\dfdf\r,s){r,s\did^\0,0) = dgd:^cf/d:,^ . (31) 

As we are dealing with freely rotating atoms, the above 
transition elements must be averaged over all directions. 
Furthermore, assuming isotropy of the atoms, we have 


fjOrJrO — X. 


(32) 


E''gl{vA,VB) = Ef{vA)EEf{vB) + 

+ ELon{T^A,YB) + ENAi'TA^YB) , (26) 

where 


where the symbol denotes the average over all direc¬ 
tions. An analogous equation holds for atom B. So, 
collecting the last two equations and substituting them 
into (1301) we obtain London’s result in its most common 
form|41|. 


Eloti — ^ ) 


, {0,0\HAB\r, s){r,s\HAB\0,0) 
Eg + Eg 


(27) 


and 


(0,0|RAB|r,s)(r,s|R^Bs|0,0) 


Ena = — ^' 

'^Or -‘-^Os 

(0,0|RABs|r,s)(r,s|RAB|0,0) 


Eg + Eg 


-E' 


Eg + Eg 


(OjOlRABsIr-, s)(r, s|i7Ass|0,0) ^ 

^ ^ ‘ ’ 


Hence, we are left with just four terms to analyze. The 
one contained in Ebou is independent of the surface and 
the only one that survives in the absence of the surface 
{Gh =0), so it accounts for the interaction between the 
atoms in the vacuum. The other three terms, expressed 
in Ena, are the key players in this paper and contain 
the non-additivity effects. Before delving deeper in this 
contribution, though, let us focus on Ebou in order to 
re-write it in a more familiar way. From equation m 
we see that 


EIab 




(dfdf)V;V,-: 


r—rA,r—rB 


Elotl 


_1_W' 


(Eg+Eg) • 


(33) 


Having re-obtained this important expression as a partic¬ 
ular case of our general expressions, we turn to the non 
additive terms in next section. 


III. THE NON-ADDITIVITY TERM 


Let us now concentrate our attention in the non addi¬ 
tive effects. Equation (1^ is the only one that depends 
simultaneously on both atoms and the surface - it reflects 
and contains the non-additivity effects that are inherent 
to the van der Waals dispersive interaction. The math¬ 
ematical treatment of this term is completely analogous 
to that given to Ebou in the last section. In such a way, 
the first two terms on the right-hand-side of ea. (E51) are 
equal and given by 


-E 

r,s 

= -E 


, (0, 0\HAB\r, s){r, s|iJAB5|0 > 0) _ 
Eg + Eg " 

, (0, 0\HABs\r, s){r, s\Hab\0, 0) _ 
Eg + Eg 

? iEg+Eg)r^^^ 


367regR^5 
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— 3{RAB)i{RAB)jGj^ {J^A,TCb) 


(34) 


IV. APPLICATIONS 


while the last term on the right-hand-side of ea. ((28ll is 


A. Two atoms and a conducting infinite plane 


-E 


, (0 , OlifASsk, s|^^ylBs|0,0) 


K + K 

jor 1 2 1 jOs |2 


1 !rH( xi2 




where we defined 

a"(ri,r2)=V,V'Gj^(r,r') 

r=ri,r'=r2 

Substituting eqs. (I34II35I) into ea. (l28l) we obtain 


(35) 


(36) 


As a first example let us consider two atoms in front an 
the influence of an infinite conducting plane as illustrated 
in figure [T] 



Ena = E^^\ + , (37) 


with 


E 


( 1 ) _ 
NA — 


Rab 

ISttEqA^ 


{va,vb)+ 


— i{RAB)i{RAB)jG^j{YA,VB) 


E 


( 2 ) 

NA 


Rab 

'del 




where 


Rab = 

r,s 



(38) 

(39) 


(40) 


Eos. (1551) and (1551) constitute the main result of this paper. 
We now proceed to evaluate the influence exerted by the 
surface on the interatomic interaction by analysing the 
ratio 


Ena 

El on 


AttR\b 

3 


TrC/'^(rA,rB) — S{RAB)i{RAB)j 


X gE{rA,rB) 


8^^R\b 

3 


E 


2 


g^{rA,rB) 


(41) 


The influence of the surface on the atomic interaction, in 
this order of perturbation theory, is a purely geometrical 
effect, since it doesn’t depend on the internal structure of 
the atoms. The formulas developed in this section show 
that in order to calculate the non-addivity effects all we 
have to do is to find Gh, which can be done by mapping 
our problem into an electrostatic one, according to the 
eq.(|ni). In the next section we will illustrate this method 
with treating some examples. 


FIG. 1: (color online) Two atoms, A and B, in the presence 
of an infinite conducting plane. A is the image of A and the 
image of B is not represented. Rab is the distance between 
the atoms and Rab is the distance between B and the image 
of A, denoted by A. 


We choose the axis as to have the conducting plane at 
z = 0 and both atoms belonging to the XZ plane. By 
eq. dH) we get 


G//(r,r') = - 


Xk\J{x — x'Y -\- [y — y'Y + {z + z'Y 


(42) 


Employing eq. dMl), the only non vanishing terms are 


Gxx 

Gyv 

Gzz 


3 sin^ e - 1 


47ri?^^ 

1 

47ri?^^ 

1 — 3 cos^ 6 
An Rab 


^ 3 sin 0 cos 9 ^ 

tJxZ - o - ZX 5 

AnRAB 


(43) 

(44) 

(45) 

(46) 


where we used Rab sin 9 — Rab sin 9. Substituting these 
expressions into eg. (1551) we arrive at 


E 


( 2 ) _ 
NA ~ 


Rab 


2An‘^£lR 


= Elou^R) ■ 


O^^AB 


(47) 


where we used eq. (155)) . This term stands for the London 
interactions between an atom and the image of the other. 
To evaluate the other term, note that Rx = sine and 
Rr = cose, while the u-component vanishes. Hence eq. 
dMl) yields 


771 ( 1 ) __ 

^NA — 


Rab 


727r%2i?37f 


(2-3 sin^ 0-3 sin^ 9) . (48) 
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We have now the complete expression to the non¬ 
additivity terms up to second order in perturbation the¬ 
ory for the two atoms-conducting plane case. The sum 
coincides precisely to the result obtained by 
Power and Thirunamachandranjs^ . 


Working directly with G, the four terms depicted in eqs. 
O, (131 and (I35l) can be naturally assembled in just 
one term 

ELon + EMA = -^Y.'^Q^{vA.rB)]^ , (51) 

7,. 7 


B. Two atoms inside a plane capacitor 


where 


This constitutes the main example of this paper. As we 
show, the non-additivity effects can are noticeable even 
at large interatomic separations as it strongly suppress 
the interatomic interaction. 


I 



FIG. 2: (color online) Two atoms, A and B inside a perfectly 
conductor plate conductor. We choose the the plane 2 = 0 
midway between the plates. Without loosing generality, atom 
A is put at (0, 0, za) and B at {xb,0, zb)- 


It is convenient to choose our coordinates in such a 
way as to have the plane 2 = 0 equidistant to the con¬ 
ducting plates, which are separated by a distance D. We 
orient the axis to have both atoms in the XZ plane as 
illustrated in figure d The Gh is obtained through the 
solution of the electrostatic problem of one charge in the 
presence of a plane capacitor, which can be done again 
employing the image method. In this case, however, we 
must deal with an infinite series of images. The poten¬ 
tial generated by this series, however, is very slowly con¬ 
vergent and it is more convenient to write it in another 
form. We follow Ref.jl^ to find that for this geometry 
the Green function may be written as 


G(r,r') 


I 

ttH 


E 


n-KZ 

cos cos 


niTz' 

D 


Ko 


nMp- P'\ 
D 



where p = xx -|- yy and p' = x'x -I- y'y and Kq is a mod¬ 
ified Bessel function of the second kind [ 3 ^. Note that 
in the last equation figures the complete Green function 
given in eq. (ED, instead of Gh- In this example it pays 
off to work directly with G, and in the end we turn to 
Gh in order to isolate the non-additivity contribution. 

The asymptotic expansion of eq. (Hi), valid in the re¬ 
gion Ip — p'l > D, is given by 


, i / 8 nz TTZ _ 

G(r, r ) = — W 1 -- cos — cos —-e 

^ ' 47rl/ p - p'\D D D 


. ( 50 ) 


Qij{rAXB) 


47r ^ 


7r2 TTz' ^\p-p'\ 

cos — cos ——e n X 
D D 



r=rA,r'=rB 


(52) 


Hence, when the atoms are separated by a distance of the 
order of D or larger, the non-addivity effects shield the 
atoms exponentially. This is a remarkable result and it 
is particularly interesting when both atoms are equidis¬ 
tant from the plates. In this case, we see immediately by 
symmetry that the force exerted by the surface in each 
atom separately vanishes. Nevertheless, the plates leave 
their footprint on the interatomic force by suppressing 
exponentially an interaction that would fall with I /R\h 
in vacuum. The results we display here for two atoms 
remain valid up to second order of perturbation when 
several atoms are present. This is so because in this or¬ 
der no term in the interaction energy of the system can 
couple more than two atoms. This can be seen by em¬ 
ploying a reasoning analogous to that we used to show 
the terms in eqs. (El and El vanish. Hence, if a gas 
is rarefied enough so that its atoms are on average sepa¬ 
rated by a distance D or more, then up to second order 
of perturbation the capacitor shields the interatomic in¬ 
teraction producing an ideal gas behavior. Even when 
this condition is not strictly satisfied our results show 
that the atoms interact only with the atoms closer than 
H, leading us to the expectation that a gas put between 
parallel conducting plates behaves more ideally than oth¬ 
erwise. Since the interatomic interaction is the main fea¬ 
ture responsible for the gas-liquid phase transition, we 
conjecture that putting a gas between conducting plates 
could lower the liquefaction temperature. 

At last, we isolate the non-additivity effects through 
equations (l38l) and (l3^ . In order to obtain Gr we must 
subtract I/|r — r'| from ea. (H51) . In figure [3] we plot the 
non-additivity effects normalized by the London energy 
for both the complete expression and the asymptotic ex¬ 
pansion (obtained using ea. dSOt for Gh), in the situation 
where both atoms are 2 = 0 (equidistant from the con¬ 
ducting planes). As expected from our previous discus¬ 
sion, for large distances the ratio goes to —I, showing 
that the non-additivity cancels out the London interac¬ 
tion between the atoms. Finally, in figure E] we do the 
same thing for the p-component of the force on atom A. 
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Ena 



^AD 


FIG. 3: (color online) Non-additive part of the interaction 
energy (normalized by the London ineraction energy) as a 
function of the distance between atom A and atom B. The 
atoms are both equidistant to the plates. Rab measured in 
units of D. The dashed bine curve employs the asymptotic 
expression given in eq. (O and the red solid curve employs 
the complete expression (f49l). 


^NA 



Ead 


FIG. 4: (color online) Non-additive part of the force (normal¬ 
ized by the force between the atoms in vacuum) as a function 
of the distance between atom A and atom B. The atoms 
are both equidistant to the plates. Rab measured in units of 
D. The dashed blue curve employs the asymptotic expression 
given in eq. dsni and the red solid curve employs the complete 
expression dm. 


C. Two atoms and a conducting grounded sphere 

In this section we analyse the non-retarded interaction 
between two atoms and a grounded perfectly conducting 
sphere. Since the application of the image method to 
the sphere is well-known, we shall use it to obtain Gh- 
Placing the origin of our coordinates at the center of the 
sphere of radius a and setting r' as the position of the 
physical charge, we have to put an image charge qi = 
— frq at position r' = Therefore, from eq. (0 

we have 


47 rr'|r —r/| dTry/r^r'^ — 2r • r'o^ -|- 

(53) 


Following the scheme outlined in the last sections, we 
now use this function to analyse the quantum dispersive 
interaction between the atoms and the sphere. Calculat¬ 
ing Qij from eq. (1551) . we obtain after some algebra 


g 


H 

ij 


3a{xfrg - xfa^){xfr\ - xfa^) 

ATT[r\r'^g — 2 va- 

a{2xfxf — 

47r[r^rg — 2 ya ■ ’ 


where xf stands for the i-th cartesian coordinate of 
r^, rA = ItaI, with analogous notation for the coordi¬ 
nates oi B. It is convenient to orient the axis in or¬ 
der to have the atom A at (0,0, ta) and the atom B at 
(0,rB sin0,rB cos6(). Employing equations (1551) we get 
the complete expression for the non-additivity terms for 
any configuration of the atoms, but for the sake of clarity 
we write explicitly only two particular cases. When the 
two atoms are aligned with the center of the sphere, we 
see from ((54l) that Q is diagonal. In that case, we have 
0 = 0 (0 = tt) when the atoms are on the same side (op¬ 
posite side) of the sphere, thus the non-additivity terms 
are given by 


B 


( 1 ) 

‘NA 


E 


( 2 ) 

NA 


Aab arATB , . 

(j'acb ± 

Aab 3a^ =F 2a'^rArB + a'^r\r% 
Idd-TT^EQ {ta^b ± a?Y 


where the upper (lower) sign refers to the 0 = n {0 = 0) 
case. Note that both terms are bigger for 0 = 0. Since, 
by symmetry, we expect the same behaviour when we 
change 0 by 27r — 0, we conclude that 0 = tt will be a 
minimum for the non-additivity interaction energy. 

In this case, in contrast with the capacitor example, 
non-additivity effects are practically unnoticeable since 
by far the dominant contribution comes by far from the 
attraction of each atom with the sphere, which con¬ 
tributes in first order of perturbation theory. However, 
since symmetry considerations ensure that the force ex¬ 
erted by the sphere on each atom separately is radial, 
non-additivity may be relevant in the components of the 
force transverse to that direction. Let us then consider 
the situation where the atom A is at the 2^-axis and Rab 
is parallel to the y direction, as illustrated in figure [5] 

In this setup, only the London interaction and non¬ 
additivity terms contribute to the j/-component of the 
force on atom B. In figure [5] we plot this component of 
the force as a function of the distance from atom B to the 
center of the sphere, keeping the vector Rab fixed. We 
see that the non-additivity parcel may indeed be com¬ 
parable to London interaction for close distances. For a 
sphere of 1/im of radius, and the atom R at a distance 
of 1 nm from the surface of the sphere, and separated by 
a distance 2 nm from atom H, the non-additivity force is 
30% from the London force between the atoms. 

To end this example we make some last remarks. First, 
remember that in our first example of two atoms in the 


















FIG. 5: Two atoms near a conducting sphere. Ras is parallel 
to the T-axis. 


Fy 

Flm 



FIG. 6: (color on-line) Non-additive part of the x-component 
of the force exerted on atom B (normalized by the London 
force) as a function of the distance between atom B and the 
center of the sphere. Ras remains always perpendicular to 
vb- The red solid curve is for two atoms separated by dA = 
0.002a and the dashed blue curve is for dA = 0.003a. 


presence of a plane, given in ea. dTfl) . could read¬ 

ily be identified as the London interaction between each 

( 2 ) 

atom and the image of the other. In this example, ifjvA 
does not have such a simple interpretation. This is re¬ 
lated to the fact that to solve the electrostatic problem of 
a dipole in the presence of a grounded conducting sphere 
we must have not only an image dipole but two point 
charges as well [ 3 ^. As a particular case of our results we 
may obtain, by substituting ea. (l53l) into ea. dT^ the dis¬ 
persive interaction between one atom and a conducting 
sphere. In doing so we arrive at the same result obtained 
in Ref. HI. 


D. Two atoms and a conducting isolated neutral 
sphere 


We need only minor modifications to tackle the case 
where the two atoms are in the presence of an isolated, 
instead of a grounded sphere. In the electrostatic case, 
the isolated sphere interacts weaklier with a charge, since 
in the grounded case the sphere is supplied by the earth 
with additional charges. Let us see what happens to the 
non-additivity effects in this quantum problem. 

Gh is obtained from the one in the last case by adding 
only one termj^. namely, 

Gh{y, r') =- , ” -- + . (57) 

47r\/r^r'^ — 2r • r'o^ -I- 47rrr' 

Substituting eq. (EZD into eg. (1361) and using our previous 
result, eq. (IMl) . we obtain 


Mxfrl - xfa^){xfr\ - xfa^) 
47r[r^rg — 2rA • -I- a'^Y’G 

aj^xfxf - Sija^) axfxf 

ATT[r\rQ — 2rA ■ 


Once more, we may evaluate the non-additivity contri¬ 
bution to the interaction energy of the system in any of 
its configurations. As before, we write explicitly only the 
case when the atoms are co-linear with the center of the 
sphere. Using eqs. (|55)) . we obtain 


Fna — 

£’0) _ _ 
^NA — 


Aab 


36tt^sIR% 


arATB 


{tatb ± 


^A^B 


A 


AB 


1447r%2 


/ aTATB T cY _ 
(tat 

2a^ 


\{rArB ±a'^Y 


{vatb ± 0^)6 


(59) 


(60) 


It can be seen in figure [7] that the non-additivity effects 
are greatly diminished in comparison with the previous 
case, where the sphere was grounded. 

Let us examine an interesting particular case. Letting 
a —)> 0, we see that the first contribution to 
is of third (fourth) order in a. Therefore, keeping only 
the lowest order contribution, we may write for the non¬ 
additivity term of two atoms plus a very small sphere the 
expression 

JfwA = - a a^ + 0(aY, (61) 

YYab' A' B 

where AT is a positive constant. In this limit we may 
identify the shrunk sphere as an atom with polarizability 
a{G) proportional to a^. Also, ta and tb become the 
respective distance from the atoms A and B to the atom 
G. Therefore, we see that our result is compatible to the 
well known Axilrod-Teller potential for three atoms, with 
the expected negative signQ. Note that such compati¬ 
bility follows from the absence of the powers a°, a and 
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FIG. 7: Ratio between the non-additive the isolated sphere’s 
and the grounded sphere’s non additive contributions for the 
interaction energy. The horizontal axis stands for the distance 
of atom A to the center of the sphere (in units of the radius of 
the sphere). The atom B is farther from the sphere, separated 
by a distance Rab = 0.002a from atom A and co-linear with 
it and the center of the sphere. 


in the above expansion, which didn’t happen for the 
grounded sphere in the last section. As it should, since 
this limit on the grounded sphere can by no means be 
thought of as an atom. We may yet present eg. (15^ in 
a more familiar form. For two-level atoms we havefdH 
Aab proportional to a{A)a{B), leaving eq. (15^ in the 
form 


Ena = — 


Ka(A)a{B)a{C) 


R?r\ 


O(a^). 


( 62 ) 


A' B 


V. FINAL REMARKS 

We have dealt with two atoms in the presence of a 
conducting surface of an arbitrary shape. In systems 
composed by three bodies it is well-known that non¬ 
additivity effects must be taken into account. In the case 
of three atoms the effects are rather small since they are 


of third order in perturbation theory while the additivity 
terms are of second order. When a conducting surface 
is present, however, the non-additivity term is also of 
second order and may be relevant to the interaction be¬ 
tween the atoms. We have obtained explicit expressions 
that allow us to evaluate analytically such an influence 
from the solution of a considerably simpler electrostatic 
problem of only one charge with the conducting surface, 
which may be treated by image method whenever pos¬ 
sible. We checked the self consistency of our results re¬ 
obtaining London formula. Besides, for the case of two 
atoms and a conducting plane we also re-obtained the re¬ 
sult displayed in literature. Then we discussed our most 
important example of two atoms inside a plate capaci¬ 
tor, where the non-additivity can not be neglected. We 
showed that the non-additivity shields one atom from the 
other, making the interaction between them to fall expo¬ 
nentially with the distance. This effect is present also 
for N atoms between two inhnite planes up to second or¬ 
der in perturbation theory. In such a way we concluded 
that a gas is closer to an ideal one between conducting 
planes, leading us to conjecture that the gas-liquid tran¬ 
sition takes place at lower temperatures inside a plate 
capacitor than inside non-conducting plates. As a last 
example we treated two atoms in a presence of a con¬ 
ducting sphere both grounded and isolated. We demon¬ 
strated that when the sphere isolated the non-additivity 
is much smaller than in the grounded case. The isolated 
case, however, has a nice particular limit, namely, the 
three atoms configuration, obtained when we let the ra¬ 
dius of the sphere go to zero. We expect that the general 
and simple nature of the results enlisted in this paper 
allows for a broader understanding of non-additivity, in 
situations where the distances involved are small enough 
to a non-retarded treatment be appropriate. 
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